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In the wide variety of axion cosmologies in which axion strings form, their radiative decay is the dominant 
mechanism for the production of axions, imposing a tight constraint on the axion mass. Here, we focus on the 
mechanism by which axions are produced in this scenario and, in particular, the key issue of the axion spectrum 
emitted by an evolving network of strings. 


1. INTRODUCTION 

The spectrum of radiation emitted by global 
cosmic strings is of crucial importance to viability 
of axionic cold dark matter scenarios. There are 
a broad variety of scenarios in which a network of 
cosmic strings are formed in the early universe, 
irrespective of whether or not in an inflationary 
context the reheat temperature of the universe is 
above the Peccei-Quinn scale /a (see ref. |p^). 
The axion string network evolves quickly to a 
scaling regime in which there are a fixed num¬ 
ber of long strings stretching across each horizon 
volume. This overall energy loss is maintained by 
the production of loops and subsequent emission 
of massless axion radiation. Around the QCD 
phase transition, ‘soft’ instanton effects break the 
residual discrete symmetry and the axions be¬ 
come massive. This process also produces domain 
walls connected to the strings; a hybrid system 
which is unstable and collapses within one Hubble 
time, if there is a unique CP conserving vacuum 

iN=l). 

The important issue when considering the con¬ 
tribution of axions to the mass density of the uni¬ 
verse is the effect of high frequency modes and the 
literature on this subject has traditionally split 
into two unequal camps. The original work [Q 
assumed that the spectrum was dominated by 
modes close to the fundamental frequency of 
the perturbation, a contention which has been 
broadly supported by numerical simulations of 
axion radiation from an oscillating string , plus 
a detailed comparison between more extensive 


numerical simulations and analytic predictions 
treating the strings as line-like objects within the 
framework of the Kalb-Ramond action ||]. Un¬ 
der these assumptions one would see a perturbed 
string oscillate, losing energy all the time, with 
the maximum amplitude of oscillation being grad¬ 
ually damped. If one assumes that the domi¬ 
nant mechanism for the production of axions is 
through the collapse of string loops then 

the prediction obtained by assuming Ua ~ 1 is 
that the axion mass ma ~ 100/reV, with possi¬ 
ble uncertainties, most of which are common to 
all calculations of the cosmological axion density; 
these are about a factor of 5 in each direction (see 
ref.@). 

This widely accepted viewpoint has been re¬ 
peatedly questioned by Sikivie and collabora¬ 
tors 1^,^,^, who maintain that the power emit¬ 
ted in the nth mode falls-off like P{n) oc n~^ 
for large n, that is, a flat spectrum per logarith¬ 
mic frequency interval. A physical consequence of 
this prediction is that a perturbed string would 
straighten itself in one oscillation, that is, it would 
be critically damped. This divergent spectrum 
leads to a prediction which systematically differs 
from the standard prediction by up to two orders 
of magnitude Wa 1/ieV with similar levels of 
uncertainty. 

These two very different predictions for the ax¬ 
ion mass can be thought of as being due to one 
universal difference between the predicted spec¬ 
tra. If one considers a spectrum of radiation 
P{n) cx n~'^, then the calculation of the axion 
density involves an integration of this spectrum 
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over the range n = ni to n = 712 , where the fre¬ 
quencies Til oc 1/t and 712 oc /a correspond to the 
long distance and short distance cut-offs, natu¬ 
rally given by the radius of curvature of the string 
(~ t) and its width (~ /a”^)- In the case where 
q > 1 then the dependence on the high frequen¬ 
cies is proportional to {n 2 /niy^^ (plus an 0(1) 
factor dependent on q) and, therefore, in all realis¬ 
tic cosmological scenarios this dependence is weak 
and effectively can be ignored. However, the sit¬ 
uation is very much different if <7 = 1. In this case 
the integral is logarithmically divergent - hence, 
it is called a ‘flat spectrum’ - and the axion den¬ 
sity is proportional to log(7ii/772) ~ 100. At the 
simplest level, it is this factor 100 which is respon¬ 
sible for the two orders of magnitude discrepancy 
in predictions for the axion mass. Hence, from the 
point of view of axion cosmology the important 
question is which is the better approximation: (i) 
a spectrum dominated by the lowest harmonics 
q > 1 or (ii) the flat spectrum q = 1. 

In the rest of this article we shall present ana¬ 
lytic and numerical arguments which clearly sug¬ 
gest that q > 1. First, we will summarize the ba¬ 
sic points of work presented in ref. which con¬ 
firms the efficacy of the Kalb-Ramond (KR) ac¬ 
tion for describing the dynamics of global strings 
by comparison to numerical simulations. In do¬ 
ing so we will respond to various criticisms which 
have been made about our methods, showing 
them to be groundless. This work supports very 
strongly the basic physical premises of the orig¬ 
inal picture for global string radiation put for¬ 
ward by Davis. Then we will report on new 
work in which we evolve a small network of global 
strings |j^. Although these simulations are the 
largest ever done (300^ points), they are still very 
small relative cosmological scales and the results 
must be correctly interpreted. Making sensible 
extrapolations to cosmological scales, they also 
confirm unequivocally the view that q > 1. 

2. RADIATION FROM OSCILLATING 
STRINGS 

The basic features of global cosmic strings can 
be described by the Goldstone model with a 
global U(l) symmetry for a complex scalar field 



Figure 1. Axion radiation from an oscillating pe¬ 
riodic string configuration in a three-dimensional 
field theory simulation. In this perpendicular 
cross-section, the string oscillates horizontally. 
Note the dominance of the n=2 quadrupole ra¬ 
diation 


d), whose action is 



- Ia ($$ - /f)" 


( 1 ) 


This can be transformed into the KR action for 
line-like strings coupled to an antisymmetric ten¬ 
sor field 

S = —qiQ j ^/—^dadr ^ J 

-2ttU j B^,d<7'^\ ( 2 ) 


by exploiting the duality between a massless 
scalar field and the antisymmetric tensor A 
detailed explanation of the derivation of this ac¬ 
tion and the notation used here can be found in 
ref. ||]. For the current discussion, all we need to 
appreciate is that single string can be represented 
within both formalisms, in the field theory case 
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by a vortex solution to the field equations and in 
the KR action by a solution to the Nambu equa¬ 
tions of motion modified to include the effects of 
radiation backreaction. If {a, r) are the coor¬ 
dinates of the string, then the equation of motion 
is 

- Xn = Ft,,, + , (3) 

where F^.,, and F^^^ are the contributions to 
the backreaction from the self and radiation 
fields respectively. The dominant divergent 
part of the self force is given by Ft,, = 
—27r/f log(A/(5)(X^ — leading to the well 

known renormalization of the mass per unit 
length 

fi{A) = Mo + 27r/2 log(A/(5). (4) 

This renormalization is critical to our argument 
based on the Kalb-Ramond action since, if we as¬ 
sume the effects of the radiation field are small, 
then the zeroth-order string solution is just that 
of the simple Nambu equations of motion and the 
effects of radiation can be thought of as being 
small perturbations. It is analogous to the mass 
renormalization of the point electron used in elec¬ 
tromagnetism, and a number of papers in the lit¬ 
erature support it . 

To proceed analytically from this point we just 
need solutions to the Nambu equations of motion 
and then the power spectrum of axion radiation 
can be computed using a formula similar to the 
quadrupole formula used for gravitational radi¬ 
ation 1^. We investigated the spectrum emit¬ 
ted by a wide range of loop and periodic infi¬ 
nite string solutions, these included the Kibble- 
Turok loops - of which the circular loop is a spe¬ 
cial case, a kinky loop, a helicoidal long string, 
a sinusoidal long string and a kinky long string. 
(Note that ‘kinks’ correspond to discontinuities 
in the tangents and velocities along the string, 
which propagate at the speed of light.) In the 
only exactly tractable case, the helicoidal string, 
the power spectrum at large n was shown to be 
exponential and various approximate arguments 
suggested that the quadrupole mode was the first 
for any exactly periodic solution. All the other 
solutions studied could be shown asymptotically 


to produce the generic q > 1 spectrum, but for 
two exceptions. First, a perfectly circular loop be¬ 
cause of its special symmetry collapses to a point 
and produces a divergent <7 = 1 spectrum. How¬ 
ever, unlike this special case, loops produced by a 
realistic network can be expected to have a spec¬ 
trum at least q > 4/3. Secondly, a periodic kinky 
long string solution can be written as the sum 
of an infinite series of sinusoidal perturbations, 
so this too has an apparently divergent spectrum 
with P{n) oc n~^ for large n. However, as we dis¬ 
cuss below, radiation backreaction rapidly damps 
out the higher frequency modes, effectively turn¬ 
ing the initial kink into a sinusoidal perturbation 
radiating into the lowest modes. 

A complementary numerical approach can be 
taken by evolving the field equations correspond¬ 
ing to the action on a discretized grid (100^ 
points). In order to compare directly with the 
analytic calculations discussed above, we concen¬ 
trated on perturbed periodic long string config¬ 
urations, in particular those with sinusoidal and 
kink perturbations. Separation of the self-field 
from that of the radiation poses a more serious 
problem in this approach than in the KR formal¬ 
ism. In fact there is no exact way of doing this, 
but it can be done approximately, the errors be¬ 
ing largest close to the string core. The method 
we use is to wait until the string becomes straight 
and then subtract out the static solution. As we 
have already noted in our discussion of the KR 
formalism, the self-field of the string obeys the 
same equation of motion as the string itself allow¬ 
ing the renormalization of the string tension (H). 
Hence, as the string moves, the self-field moves 
with it, maintaining approximately the same form 
as the static solution. Obviously, Lorentz con¬ 
tractions of the string may distort the string, but 
this will only become important at large ampli¬ 
tude when string velocities u « c. 

Perhaps the most compelling evidence for the 
validity of the subtraction comes from the results 
themselves since they are so clean. As shown 
in Fig. 1^ the radiation field for a sinusoidal per¬ 
turbation comprizes of almost perfect quadrupole 
lobes characteristic of the n=2 mode, which can 
be seen to propagate outwards from the centre 
of the string as the system evolves. It has been 
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suggested that this travelling wave phenomena 
could be an artefact of the subtraction, but this is 
simply not credible! Further confirmation comes 
from spectral analysis of the radiation field. We 
find that the spectrum is dominated by the n=2 
mode with a rapid fall-off into the higher modes; 
any problems with the subtraction would surely 
create high frequency noise. Further discussion of 
this issue will be included in ref. . 

If one sets up an initial kink solution, the spec¬ 
trum of radiation initially contains more power 
in higher frequency modes, as suggested by the 
analytic calculations. But as the solution evolves 
higher frequency radiation appears to be damped 
out much faster than the lower frequencies, close 
to the characteristic frequency of the perturba¬ 
tion. Soon the spectrum of radiation becomes 
indistinguishable from that of the sinusoidal per¬ 
turbation. The initially sharp profile of the kink 
becomes visibly rounded, so that after around 5 
oscillations the solution looks very much like a 
sinusoidal perturbation. This strongly suggests 
that high frequency radiation is damped much 
more quickly by radiation backreaction, generi- 
cally leaving behind the dominant lower harmon¬ 
ics [§. 

3. THE COLLAPSE OF A NETWORK 
OF LOOPS 

In order to go further than the single string sim¬ 
ulations discussed in the previous section, we have 
evolved a network of strings in a 300^ box with 
periodic boundary conditions. The intention here 
is not to investigate the scaling properties of the 
network, but rather to analyze the spectrum of 
radiation emitted by more general string config¬ 
urations. The initial conditions were created by 
laying down a random noise configuration, before 
initially evolving it in a dissipative regime (re¬ 
moving the radiation background). This process 
simulates some aspects of a cosmological phase 
transition. However, for our purposes the pre¬ 
cise details of the initial distribution of strings 
are not important, the only salient point being 
that a simple random network of loops (due to 
the periodicity) can be created which in many 
ways mimics that formed in the early universe. 



Figure 2. The results of 3D spectral analysis 
of the argument of the complex scalar field in 
numerical simulations. We plot \Ek\ = k‘^\9k\‘^ 
against k, which is the square root of the spec¬ 
trum of the axion energy density. On the right 
is that from the collapse of loops, included are 
lines representing k~'^ (dotted line) and 
(dashed line), and on the left from the collapse 
of strings connected to domain walls, with lines 
for k~^ (dotted line) and (dashed line). 


As the network evolves relativistically some of 
the loops collapse and emit massless axion radia¬ 
tion, while others interact with each other creat¬ 
ing larger and smaller loops. This process con¬ 
tinues until no more interactions are possible, 
with all loops collapsing into axions. The spec¬ 
trum of loop sizes creating the axions has a fairly 
wide random distribution, being a function of the 
initial conditions and the subsequent dynamics. 
Once all the strings have collapsed all that is left 
is axion radiation, so there is now no self-field as¬ 
sociated with the strings and hence we do not 
have to worry about any subtraction ‘ambigu¬ 
ities’ before performing a 3D spectral analysis. 
We have analyzed the argument of the field $, 
which corresponds to the massless axion field, ig- 
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noring the massive radiation emitted by strings 
on these small lengthscales. This is plotted as 
\Ek\ = in Fig. ||(a), which is the square 

root of the power spectrum of the axion contri¬ 
bution to the energy density. At small k - that 
is, large wavelengths - the spectrum appears to 
be dominated by the initial correlation length of 
the strings, that is, wavelengths corresponding to 
the most likely loops in the simulation (well above 
the string width). However, the focus of our at¬ 
tention is on the distribution of high frequency 
modes. For k > 10, we see that the spectrum of 
axions present can be approximated by a fairly 
steep power law |iffc| oc k~‘^. Despite the consid¬ 
erable uncertainties in the exponent (which we 
shall discuss elsewhere [^), these results are un¬ 
doubtedly inconsistent with the flat q = 1 spec¬ 
trum which should fall-off here at the much slower 
rate 

We have also performed simulations starting 
with the same initial conditions, but with the ax¬ 
ion mass switched on from the start. In the orig¬ 
inal simulations the action was rescaled so that 
the problem is dimensionless and /a = 1 in these 
units. When the axion mass is included the value 
we use is ma = 0.5, which is very large compared 
with the cosmologically interesting value. But for 
the mass to have an effect in these simulations on 
a reasonable timescale one requires rria = 0{fa). 
This should be kept in mind when interpreting 
these preliminary results. 

We find that the network collapses much more 
quickly once the mass is included since domain 
walls are formed connected to the strings. This 
hybrid network is unstable topologically, and the 
relativistic dynamics of the strings slices up the 
network creating small pieces of domain wall 
bounded by strings, which subsequently decay 
into axions. Once all the pieces have collapsed we 
analyzed the spectrum of axions left in the box. 
As one can see from Fig. ^(b) the spectrum can 
be approximated by a power law, this time some¬ 
thing close to the which is again definitely 
not consistent with the flat spectrum. Therefore, 
we have shown that in the QCD regime where 
the axions mass ‘switches on’, the spectrum of 
axions emitted may have a slightly harder spec¬ 
trum with more power in high frequency modes. 


albeit for small networks with an artifically high 
mass. However, there appears to be no evidence 
whatever for anything even close to a flat spec¬ 
trum. 

We note, however, that it would be legitimate 
to criticise these simulations on two important 
counts. First, the size of the simulations is very 
small compared to the cosmological scales we are 
trying to model. The renomalized mass per unit 
length is ^(A) « 3^0) as compared to the realis¬ 
tic value of /i(A) « 100/io, making the radiative 
decay of the strings much more rapid than one 
would expect for a cosmological network. How¬ 
ever, such a point could only have some verac¬ 
ity if we were trying to defend a spectrum which 
was close to flat, since the damping in the simu¬ 
lations is much closer to ‘critical’. Secondly, the 
mass of the axion used in the domain wall simula¬ 
tions is far higher than that for the early universe. 
However, the effect of an unphysically large mass 
should be to speed up the annihilation process, 
making higher frequency modes more prevalent 
than in a cosmological setting. We shall return 
to these issues at greater length in a subsequent 
publication |ll| ]. 

4. DISCUSSION 

In this paper we have attempted to elucidate 
the reason why the technical issue of the power 
spectrum of radiation emitted by global strings 
is important for axion cosmology and hence for 
the large-scale axion dark matter searches. At 
a first approximation for an order of magnitude 
estimate of the axion mass, all one has to do 
is decide between the two qualitatively different 
spectra with exponents g > 1 and g = 1. We 
have summarized a number of ways, both ana¬ 
lytic and numerical, by which we can gain insight 
into what the true spectrum of axion string radi¬ 
ation actually is. In each case, these approaches 
demonstrate a g > I spectrum dominated by the 
lowest harmonics, one which is far removed from 
the flat g = I spectrum required for a low axion 
mass rria ~ I/xeV. While each of these approaches 
has caveats, the consistency of the picture they 
present is compelling! Having established the un¬ 
derlying axion spectrum on a firmer foundation, 
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the time has come to unite our efforts in mak¬ 
ing accurate predictions for the mass of a dark 
matter axion by addressing the real quantitative 
issues that remain. 
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